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Abstract 

We will show that the ADM 3+1 evolution equations, for a zero shift 
vector, arise naturally from the equations for the second variation of 
arc-length. 

1 Introduction 

Discussions of the dynamics of general relativity often begin with the ADM 
3+1 evolution equations [T]. These equations describe the second time deriva- 
tives of the spatial metric in terms of other quantities such as the lapse func- 
tion and the extrinsic and Riemann curvatures. If by some means we happen 
to have a local solution (in time) of those equations then we could imagine 
computing the arc length for short geodesic segments. What then would 
we get for the value of the second time-derivative of that arc-length? This 
question has been discussed many times elsewhere |2l [3] but under a differ- 
ent name - the second variation of arc- length. Clearly the second variation 
of arc- length and the ADM 3+1 evolution equations must be related. The 
purpose of this paper is to establish that relationship. The result is not un- 
expected - the equations for the second variation of arc-length can be used 
to recover the standard ADM 3+1 evolution equations with zero shift vector. 

There is value in this presentation beyond the purely pedagogical - the re- 
sults presented here provide strong theoretical support for an approach to 
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numerical relativity being developed by the author [H El El [7j . This method is 
known as smooth lattice relativity and is closely related to the Regge calculus 
[HI El [in]- Both methods use a lattice to describe the metric but they differ 
most notably in the way they treat the curvatures. In the Regge calculus 
the metric is piecewise ffat with the result that the curvatures are distribu- 
tions on the 2-dimensional subspaces known as bones (or hinges) while on a 
smooth lattice we allow the metric to vary smoothly in the neighbourhood of 
any vertex. This allows all the usual tools of differential geometry to be ap- 
plied to the smooth lattice. In particular we can easily compute the Riemann 
and extrinsic curvatures in terms of the geodesic arc-lengths of the lattice 
and thus, using equations ( 3.3[3.4 ) or (2.5), the second time derivatives of 
the leg-lengths. This makes the study of dynamics on a smooth lattice quite 
simple in principle (though as with any numerical method in general relativ- 
ity the practical aspects are far from trivial). Attempts have been made to 
adapt the ADM 3+1 equations to the Regge calculus [HI [I2] but progress 
has been slow. A much more promising scheme, for the Regge calculus, is 
due to Sorkin [13] with later development by Barrett et al. [H] and Gentle 
and Miller 



2 First and second variations 

Discussions on the first and second variations normally arise when asking 
questions about geodesies such as: Is the geodesic that joins two points 
unique? Is it the shortest geodesic? How far can the geodesic be extended 
before it fails to be the shortest geodesic? The mathematical theory that 
answers these questions is very elegant and has previously found its way into 
general relativity as a tool in studying the global properties of spacetime [I6] . 
Hawking and Penrose [T71 [18] made extensive use of the first and second vari- 
ations of non-spacelike geodesies in their singularity theorems. In contrast, 
our interest in the first and second variations is that they provide a natural 
setting in which to ask different questions of geodesies: How can the first 
and second time derivatives (of the arc-length) be computed? And how are 
they related to the curvature tensors? As already noted in the introduction 
these questions will lead to the standard ADM 3+1 evolution equations with 
a zero shift vector. But first we need to introduce some basic notation and 
to make clear the class of curves we will be working with. 

Choose a point i and a small neighbourhood of i in which the spacetime is 
non-singular. All of the curves we are about to construct will have a finite 
length and will lie totally within this neighbourhood. Through i construct a 
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timelike curve Ci with affine parameter rj. From Ci we can construct a nearby 
curve Cj by dragging Ci sideways a short distance (i.e. drag Ci along a short 
spacelike vector field defined on Cj). We now have two nearby timelike curves 
Ci and Cj (the point j on Cj can be easily identified - it has the same r] value 
as i has on Cj). We will assume that the two curves Cj and Cj are sufficiently 
close that, for any given t], we can construct a unique geodesic that joins the 
two curves. 

Consider now the family of geodesies generated by allowing rj to vary. This 
family of geodesies (actually, segments of geodesies) will cover a 2-dimensional 
subspace (like a taut ribbon) which we will denote by S. We will introduce 
coordinates on in a rather obvious way. Consider a point P on S. There 
will be exactly one space like geodesic (of S, by assumption) that passes 
through P. The point P will be located some fraction. A, along the geodesic 
from Cj to Cj. Thus the points on Cj will have A = while those points on 
Cj will have A = 1. We will take the other coordinate for P to be the value 
of 77 that identifies this geodesic from all others (in S). The coordinates for 
P are then taken to be (A, 77). This situation is displayed in figure ([T|. 

Consider now a global coordinate system, x^, for the spacetime. Then S can 
also be described by functions of the form x^{\,rj). We now define a pair of 
vectors r]^ and A'^ by 

r?'^ = — = ^ (2 1) 

and a pair of unit vectors ra^ and by 

= I^V' , = ^A'^ (2.2) 

where and M are scalar functions that ensure that the vectors are indeed 
unit vectors. Clearly the vector r^^ is tangent to the A = constant curves 
while A^ is tangent to the rj = constant curves (and both vectors will, in 
general, be neither unit nor orthogonal, despite appearances in figure ([T])). 
It is rather easy to show that M = ds/dX = Lij where s is the proper 
distance along the geodesic and Lij is the length of that geodesic. Recall 
that ds/dX is constant along a geodesic while Lij = (ds/dX) dX and thus 
Lij = ds/dX. Next, using the requirement that be a unit vector leads 
immediately to M = ds/dX = Lij as claimed. Later, when we specialise to 
the ADM 3+1 formulation in section |3] we shall see that is the usual lapse 
function associated with the time coordinate r]. 

We can now state clearly the equations for the arc-length and their variations. 
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Arc-length 



U,= I [9,.^^) dX (2.3) 



First variation 

dLi 



drj 

Second variation 

d^L, 



rj 

[m^ri^y. = I m^m^ri^,^ ds (2.4) 



drj 



+ f (r?^;ar?^.m"m" - (m^m^T^^;,)') ds (2.5) 



For ease of reference we have included a proof of the above equations in the 
appendices. See also [21 [3] for more details. 



3 The ADM evolution equations. Pt. 1 

Consider a typical Cauchy surface S and suppose that the pair of time like 
curves Ci and Cj intersect E at the points i and j respectively. At i we 
have two vectors n^, the unit normal to S and m^, the unit tangent to the 
geodesic that connects i to j. If we construct a third unit vector m'^ as a 
linear combination of and m^, 

m'^ = cosh 6 + sinh 6 



we can, by careful choice of the boost angle 6, ensure that m''^ is tangent to 
S. That is we require 6 such that = n^m''^. This arrangement is shown 
in figure ([2]). In what follows we will be looking at the behaviour of various 
expressions in the case where Lij is small. So our present task is to ask : How 
does 6 vary with Ljj? The first observation is trivial : ^ as L^j 0. 
Now from = n^m'^ we have 

sinh 9 = n^m^ cosh 9 
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and thus across the leg we have 



[sinh Oy. 



[ri^m^ cosh ( 



at] ^ ' 

If we choose Ljj to be sufficiently small then we can be sure that the geodesic 
(that joins i to j) intersects S only at i and j (see figure ([s])). From this 
constraint we observe that Qi and dj must be of opposite signs and thus 



I sinh I + I sinh 6',- 



drj 



Thus each term on the left must be of order O [dLij/dt]), that is 

'dLi 



9 = 



drj 



as Lj 







3.1 The first variation 

Our aim in this section is to recast the expressions for the first and second 
variations in terms of the familiar ADM data, the lapse, shift and extrinsic 
curvatures. 

The extrinsic curvature, K^^, can be defined in a number of ways (see [T]), 
such as 

NK^, = -Nn^,, - ±(iV,^)n, 
where _L is the projection operator = h^u = + n^n^). Then 

and thus 

= (m^A^,^) sinh^ - NK^.^m'^m'' 



This can now be substituted into the integral for the first variation (2.4) 



dLij 
drj 



m^m!'7]^.yds 
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Recall that we are dealing with short geodesic segments. Thus we can use 
any of a number of methods to estimate the integral. To be specific, we will 
chose a mid point rule (see [12]) which leads to 

^ = {m^N^^ sinh^) - {NK^^m^'rrf) U-j + O (L^) 

where each term is evaluated at the mid-point of the geodesic. But since 
9 = [dL/dr]) we see that the first term is of order O (L^) and thus 

^ = - {NK,,m^mn L,, + O (L^) (3.1) 

Notice that m'^ is the unit tangent vector at the mid-point of the geodesic 
that joins i to j and thus we have 

m^'Lij = x';-xf + (3.2) 



So, if we put Axf- = — we can rewrite (3.1 ) as 



dL^ 



= -2iVir^,A4Axr,. + O [L") (3.3) 



3.2 The second variation 

Once again we use the basic definition of the extrinsic curvature to express 
the terms appearing in the second variation in an ADM form. We will do 
the calculations by splitting our previous expression for the second variation 



(2.5) into the following pieces 

dj^LjA 



J1 + J2 + J3 + J4 



d'T]'^ 

Jl = [^°;M^^"^a] • ds J^ = - R^a^fsm^'m''r]''r]^ ds 
•^3= Vfi;aV^;um"m'' ds J4 = - / [m^nfT]^.^) ds 

J i J i 

3.2.1 The second term 

We start with this term as it requires very little work. We simply substitute 
rj^^ = Nn^ and approximate the integral via a mid-point rule leading to 

J2 = - r R^,a^fsm^'m''r]°7]^ ds = -N'^R^^^pm''m!'n''n'^ Uj + O (L^) 

J i 
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3.2.2 The fourth term 

Here we use m^m^rj^.y = —NK^yin^m^ + O (L) (the error term arises from 
the n^m^ = sinh6' = O {L) terms). Thus we are lead to 

J4 = - r im^'m''r]f,,^f rfs = - T {-NK^^m^'m'' + O {L)f ds 

J i J i 



where we have used (3.1) in the second last line. 



The remaining terms are not so easily dealt with. 

3.2.3 The third term 

For the third term the details follows 

J i 

3 

^ {- {N^^m'^ f + N^n^.^n^.pm''m'^) ds 

' (- (iV,,m")' + (±(iV,^)n„ + K^^) (±(iV'^)n^ + a) m^m^) ds 
= r (- (A^,„m")^ + N^K^^K^'prn^m^ + O (L)) ds 

J i 

The error term O (L) in the last line arises from terms of the form n^m^ = 
sinh6' = O {L). Now we use the mid-point rule, once again, to obtain 

J3 = f Vt^^aV^^urnVds = - (iV,,m")' Li, + N^K^^K^ ^m'^m^ + O [L^) 

J i 
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3.2.4 The first term 



Finally, we turn to the first term [rj^-^rj'^mall- Using the same substitutions 
as we have used before and also using Nn^N^n = dN/drj we obtain 



Ji = [v'^;n.V^rna\ 



1 dN 



We choose to write this result as a sum of two terms each of the form [•••]] 
so that we can deal with each term separately. In the first term we have 
{l/N)dN/drj which varies slowly over the short geodesic and thus may be 
taken as a constant (plus a an error term of order O (L)), thus we have 



1 dN 
N dr) 



1 dN 



N dr] 
1 dNdL 



[<m,]^- + [r]^m,i O (L) 



N dr] drj 

For the second term we use a Taylor series expansion 

d 



[NN^f^rri' 



tj.]j 

i 



ds 



= {NN^m^).^^m''Lij + 0{L^) 

= {N^m'^f Lij + NNaprn'^m^Lij + O {L") 

The appearance of the term N.^ap is encouraging - it reminds us of the similar 
term in the ADM equations. We can improve on this situation. Notice that 
m'!" ^m^' + O (L) and thus 

Napm^m'^ = N,^pm"'m!^ + O (L) 

However, m'^ is tangent to E thus we also have 

Napm^'m^ = N\^pm"'m''^ + O (L) = N\^pm°'m^ + O (L) 

where the vertical stroke denotes covariant differentiation with respect to the 
3-metric intrinsic to S. 

Combining these two results we obtain our final estimate for the first term 
in the second variation 



1 dN dLi 



N dr] dr] 



+ {N^^m^f Lij + NN\^pm''ml^Lij + O {L") 



8 



Now we can reassemble the pieces. The result is 



d^Ljj _ 1 dN dLij 1 / dLij \\ ^ . 



dr]"^ N drj drj Lij \ drj 

+ NN\^pm'^m^Lij - N^Rf,auf3m''m''n'^n^Lij + O [L^) 

We are almost finished, we just need to do a little bit of tidying up. We 
multiply both side by Lij/N and noting that 

df]'^ \ df] ) ^'^ drf 

d ( 1 dL\^ _ 1 dNdL\ ^ 1 d^L\- 



drj \N drj J iV^ drj drj N drf 
we can rewrite the above equation as 

I, (^S) - 2iV..A.SA4 (3.4) 
+ 2JV (AV^A'"^ - R^^^^nV) AigAig + O (L') 

where we have also used Ax^- = m^Lij + O (L^)- 

For completeness, we repeat here the result we previously obtained for the 
first time derivative. 



4 The ADM evolution equations. Pt. 2 

This completes the first stage of the construction. We have successfully ex- 
pressed the first and second variations in terms of the extrinsic and Riemann 
curvatures. Our second and final stage will, among other things, introduce 
the metric tensor as a replacement for the geodesic arc-lengths. As we shall 
soon see, this is not a difficult task. The most notable change is not in the 
symbols, from Lf^- to g^u, but in the structure of the equations. We will be 
re-working an equation defined over a geodesic segment into an new equation 
defined at a point. 

Consider a typical geodesic segment with end-points i and j. The time like 
worldlines Q and Cj generated by i and j are, by assumption, orthogonal 
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to the Cauchy surfaces. Thus we can use these curves to propagate the 
spatial coordinates of each Cauchy surface forward in time. This means that 
the spatial coordinates of any point on Ci are constant along Cj and thus 
= dAxfj/dr]. 

We now introduce the metric by estimating Lij using a mid-point rule for 

/ ds, 



'dX'dX 



1/2 



MdX 



dx^ dx" 
'dX'dX 



1/2 



M + O (L^) 



where each term on the right hand side is evaluated at the mid point of the 
geodesic. But we have previously shown ( 2.1|2.2 ) and (3.2) that dx^/dX = 



m^Lij = Axfj + O (L^)- We can use this to estimate L^j as 

= g.^Ax^.Ax':^ + O {L') 

We can go one step further by noting that g^^, = h^^, — n^riy and n^Ax^- 
Lij sinh^ = O (L^) and thus to leading order in L we have 



which, when substituted into the above, leads to 



(4.1) 



It is now just a short step to the finish line. First substitute (3.3) into (3.4) 



and then (4.1 ) into (3.3) and finally take the Axfj terms out through the time 
derivatives. Then notice that the Ax^- are arbitrary and that the coefficients 
of Ax'^^jAx^j are symmetric in fiu and purely spatial. This allows us to cancel 
the Ax^ from both sides of the equations after which we can safely let L ^ 
(the details of this series of substitutions and eliminations are excluded as 
they follow very standard lines). As expected the final result is nothing other 
than the familiar ADM evolution equations with a zero shift vector 



dh 



flU 



drj 



-2NK, 



dK, 



df] 
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A The first variation 



We know that the mixed partial derivatives of x^{X, rj) must commute, thus 
we must have 

and for a symmetric connection (which we are using) we also have 

which can be re-expressed, terms of the unit vectors and m^, as 

{Nn").^ {Mm") = {Mm").^ {Nrf) (A.l) 

Finally, as the vector is the unit tangent to an = constant geodesic, we 
have 

= mi'.^wy 

and 

d'^x^ dx" dx^ 
dX" "^aA dX 

We will use the above equations frequently in the following discussions. 

Here we consider the geodesic arc-length and its first time derivative, 

dx^ dx'^ \ 







ds 


Lij 


-f 

Jo 


dX 


dLij 


d 


f 


dr] 


dr] . 
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dX^^^Jo a^^l^-aA aXj 

Note that the path x'^{X,r]) in each of these integrals is a geodesic and that r] 
is constant along the geodesic. The second integral in the last equation above 
can be readily evaluated using standard techniques (expand the 77 derivative, 
swap orders of mixed derivatives, integrate by parts and impose the geodesic 
equation). The result is 



dLij 1 
dr] Li, 



dx^ dx" 
^''"'dX'd^ 



3 

= [m^r]^]l 



where we have taken the small liberty of replacing the limits and 1 with the 
more suggestive labels i and j. This is an elegant result - it shows that for 
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a geodesic segment, dLij/dr] can be computed from data defined only at the 
end points of the geodesic. This may seem simple but it hides a significant 
complexity - the data involved can only be found by solving a two-point 
boundary value problem. 

Despite this compact and elegant form for the first time derivative, we will 
now develop an alternative integral expression that happens to better suited 
to our later calculations of the second time derivative. Consider for the 
moment the quantity Q defined by 

Q = [ m^m!' {Nn^).^ ds 

J i 

with the integration path being, as expected, an r/ = constant geodesic. We 
will now show that Q = dLij/dr]. We begin by writing ds = {ds/dX)dX = 



MdX and using the commutation relation (A.l) to obtain 



Q= r mX{Mm^).^ dX 

J i 

Now expand the covariant derivative and use 1 = m^m^ and = m^m'^;i, to 
obtain 



dr] \dX J dr] dX dr] 



Thus we have shown that 



- - ^m^ir]X = I m^m'^r]^,, ds \2A\ 



dr] 

Our challenge now is to compute the second time derivative. This proceeds 
in a manner similar to the above calculation though it is a tad lengthy. 
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B The second variation 



To compute the second derivative we need only apply d/dr] to (2.4). This 
leads to 

d^L,, d r „ „ . r d 



drf drj ' drj 



■m^m''ri^.^ ds = j — {m^rrfr]'^ .^M) dX 
im,r]^,,Xn.^V'' dX 



{m^.aV^.A'' + rn^^^;^;aA'' + m^7]^'.^Xr^^) r^" dX 

We will apply various manipulations to the three main parts of this integral 
and we will make extensive use of the geodesic equations, = m^-^m'^, the 
commutation relations, X^.^rj^ = rj^-^yX^ and the observations that is a 
unit vector along the geodesic. 

We start by splitting the above integral into three pieces 

h = J' m^,^r]^.,X''7]" dX 
j 



l2= / mX;u;aX''v'"dX 



m^r7^;,A^„r7" dX 
which we will now attempt to simplify. 
Integral Ji 

Put X'^ = TjfM and m^.^aM = {m^M).^ci — ^iiM-^a and then use the commu- 
tation rule on X^-aV°' to obtain 

h = f r]^.^r]^'.ym''m''M dX - T mf,M.^7]''r]''.ym'' dX 

J i J i 

Consider the second integral in this pair and denote it by J4. Since is a 
unit vector we can slide a factor of nigmP inside M.^, like this 

(mem^M) Vj'^rj^.i.nf dX 
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The term m^iam^ is zero since is a unit vector while the remaining term 
is ripe for a commutation operation. This leads to 



So our final expression for /i is 

J i 



Integral I3 

We step out of sequence here because one term arises in this computation 
that will be useful when we tackle the second integral l2- 

This integral is shghtly easier to work with then the first integral and it will 
give rise to the Riemann tensor. The main device used here is to swap the 
order of the second partial derivatives on rj'^-u-a balanced by the addition of 
the Riemann tensor. Thus we have 

h = mX.,,.am''r]"M dX 

' + R^pauVn m'^M dX 

= Is - R^o.upm^'m''ri''i^^M dX 
where we have introduced a fifth integral, 

h = f m^r}^',^.,m''r}''M dX 




Integral I2 + h 

As we shall soon sec, the integrand for I2 + can be combined to form a 
total derivative and thus the integration is trivial. We start by forming the 
sum I2 and 
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where ds = MdX. By careful inspection of the integrand, while noting the 
geodesic conditions, = m^-yUi^ ^ it is not hard to see that the integrand can 
also be written as {'r]'^ -^rj^Tna) .^m'^ . Thus we have 



l2 + h = J' (77";^r/^m,).^m^ds 



Our job is done, all of the integrals have been evaluated as far as possible 
all that remains is to combine the above results. This leads to 

d'^L ■■ ■ P 



This last integral can be simplified slightly by introducing 



which leads to 
d^L 



dr]"^ 
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Figure 1: This figure displays the 2-dimensional surface S constructed from 
the pair of time Uke worldhnes Ci and Cj. The curve connecting i to j is a 
spacehke geodesic with length Lij. Along these geodesies r] = constant. Note that 
the tangent vectors and m'^ are unit vectors but they need not be mutually 
orthogonal. 
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Figure 2: In this figure the lower (straight) curve is the geodesic that joins i to j. 
The upper curve (which is not shown in figure ([T])) arises from the intersection of 
the Cauchy surface with the 2-dimensional surface S. The unit vectors n'^ and m''' 
are orthogonal. Note that, in general, r] is not constant on each Cauchy surface. 




Figure 3: This is a situation that we explicitly exclude. In this case the points 
i and j are so far apart that the geodesic intersects the Cauchy surface at points 
other than i and j. In this case 9i and 9j have the same signs, contrary to the 
assumptions made in the text. 
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